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■ We construct a class of projected 50(5) models where the Gutzwiller constraint of no-double- 
^\ ' occupancy is implemented exactly. We introduce the concept of projected 50(5) symmetry where all 
0^ , static correlation functions are exactly SO(5) symmetric and discuss the signature of the projected 

SO(5) symmetry in dynamical correlation functions. We show that this class of projected SO(5) 
^ | models can give a realistic description of the global phase diagram of the high T c superconductors 

■ and account for many of their physical properties. 

; 

I. INTRODUCTION 

<n : 

Recently, a unified theory of antiferromagnetism (AF) and superconductivity (SC) has been proposed for the high 
, T c cuprates 0. This theory is based on the 50(5) symmetry between AF and SC, and offers a unified description 
I ' of the global phase diagram for this class of materials. While the theory was originally proposed as a effective field 
, theory description, it was soon realized that the 50(5) symmetry could be implemented exactly at a microscopic level 
and it can also be checked numerically in common strongly correlated models such as the t — J model 
While the phase diagram ppO|-p^| and collective excitations [ p~3|Jl4| in the SC state derived from these 50(5) models 
bear strong resemblance with the high T c cuprates, and a number of novel experimental predictions have been made 
Ji5| pj[ , the Mott insulating behavior at half-filling is a puzzling aspect which challenges the fundamental validity of 
the 50(5) models ]20|-|23[] . To be more precise, the exact 50(5) symmetry requires collective charge two excitation 
at half-filling to have the same mass as the collective spin wave excitations. This condition is clearly violated in a 
Mott insulating system where all charge excitations measured with respect to a particle hole symmetric point have 
a large energy gap of few eV, while the spin wave excitations are massless. In the original 50(5) proposal, it was 
pointed out that this situation is analogous to a easy-axis antifcrromagnet in a external uniform field, and a 50(5) 
symmetry breaking term at half-filling was introduced in order to describe this asymmetric behavior between spin 
and charge. The chemical potential also introduces a 50(5) symmetry breaking term, however, it was shown that 
these two terms could compensate each other so that the static potential governing the SO (5) superspin could 
1 still be 50(5) symmetric. 

Since the asymmetry between the charge and spin excitations at half-filling is of the order of the Coulomb energy 
scale U, the 50(5) symmetry breaking terms must also be of that order. Since there are various types of symmetry 
breaking terms, one might hope that their effects could partially cancel each other to arrive at a qualitatively correct 
picture. However, this type of cancellation is very delicate, and approximate calculations could easily lead to erroneous 
conclusions. In particular, one is interested in which physical properties could exhibit 50(5) symmetric properties in 
the limit when the Coulomb gap is taken to infinity. For example, one could ask the following questions: 

1) One of the hallmarks of the 50(5) symmetry is not only the degeneracy between the AF and SC states at a 
given chemical potential, but the approximate degeneracy among all mix-states interpolating between AF and SC, 
i.e. the independence of the ground state energy on the superspin angle. What is the potential barrier separating the 
AF and SC states at their degeneracy point in the limit U — > oo? If there is a large energy barrier in this limit, one 
would argue that the concept of 50(5) symmetry is not a useful one, at least not for quantitative calculations. On 
the other hand, if the potential barrier is finite and small in the U — > oo limit, the concept of a approximate 50(5) 
symmetry would be a useful one. 

2) Exact 50(5) symmetry predicts four massless collective modes. In the half-filled AF state, besides the two 
conventional massless spin wave modes, the exact 50(5) symmetry predicts a massless doublet of ir^ modes, with 
charge ±2. However, a Mott insulator has a large gap to all charge excitations. Therefore, it is clear that one of the 
tt^ has to be projected out of the spectrum in the limit U — > oo, say the n + mode carrying charge +2. What happens 



1 



to the rest of the Goldstone modes, the 7r~ mode carrying charge —2 and the 7r Q triplet mode of the SC state? In the 
U — > oo limit, can they all be simultaneously massless at the transition point between AF and SC? Since the pure 
SC state can only be reached with a finite doping concentration, is it possible that the Gutzwiller projection does not 
affect the 7r Q triplet mode of the SC state? 

In order to address these questions, it is desirable to construct a low energy effective theory without any parameters 
of the order of the Coulomb scale U. In this work, we construct a class of projected SO (5) models which treat the 
Gutzwiller constraint exactly and locally on every site. We use this model to answer the physical questions posed 
above and show that the answers are affirmative. In the U — > oo limit, when the Gutzwiller constraint is implemented 
exactly, the ground state energy can still be SO (5) symmetric and independent of the superspin direction. After 
projecting out the 7r + mode, all other Goldstone modes remain massless at the symmetric point. The dispersion 
relation of the collective modes bear unique signature of the projected SO (5) symmetry. Furthermore, the 7r Q triplet 
modes of the pure SC states are unaffected by the Gutzwiller projection. These properties define the concept of a 
projected 5*0(5) symmetry (p50(5)), whose properties and consequences we shall explore in this paper. 

The fundamental quantity in the 50(5) theory is the locally defined five component superspin vector n a {x) — 
(m, ri2, 7t3, 714,715) describing the local AF and SC order parameters respectively. In the nonlinear a model formalism, 
these are treated as mutually commuting coordinates and their dynamics is given by their conjugate momenta p a (x) = 
(pi,P2,P3,P4,P5)- The charge operator is the angular momentum in the n\ — 71 5 plane: 

Q{x) = Lm = nip 5 - n 5 pi (1.1) 

Implementing the Gutzwiller constraint corresponds to requiring 

Q(x) < (1.2) 



for every local 5*0(5) rotor. From equations (1.1) and fll.2j ) and subsequent discussions, we shall see that the Gutzwiller 
projection in the 50(5) formalism corresponds to going from a fully symmetric 50(5) rotor model to a chiral 50(5) 
rotor model, where both the static potential of the individual rotors and the coupling between the rotors are still 
50(5) symmetric, but the rotors are constrained to rotate only in one sense in the n\ — 715 plane, consistent with 



(1.2). This observation reveals a deep connection between the Gutzwiller projection and the Lowest-Landau-Level 



(LLL) projection in the fractional quantum Hall effect [E4| . To be more precise, the Gutzwiller projection represented 



by equations (1.1) and (1.2) is analogous to the LLL projection, where all states in the LLL have a definite sign 
of angular momentum. The LLL projection can be analytically implemented by separating the cyclotron degrees of 
freedom from the guiding center degrees of freedom, which amounts to changing the commuting property between 
the X and Y coordinates to a canonically conjugate commutation relation: 

[X,Y]=il (1.3) 

where lo is the Landau length. Exploiting this analogy, we find that the original 50(5) model can be fully Gutzwiller 
projected without changing its form, if one imposes the simple quantization condition between the superconducting 
components of the superspin vector: 

[ni,n B ]=</2 (1.4) 

In the symmetric 50(5) model, the wave function of the 50(5) rotors are functions of the local coordinates n\ and 
715, while the projected 50(5) model only depends on their holomorphic combination z — n\ — in^ and is independent 
of their anti- holomorphic combination z = 71 1 + in.5. This way, we arrive at a natural projection of the 50(5) model 
where the local Gutzwiller constraint is taken into account exactly, and the resulting model is free of the large Coulomb 
U parameter. Because the functional form of the symmetric 50(5) model remain the same and only the quantization 
condition is modified upon projection, many important properties associated with the 50(5) symmetry remain. The 
central hypothesis of the 50(5) theory is that this projected model is quantitatively accurate in describing both the 
static and dynamic properties of the high T c cuprates, and we shall compare the properties of this model with the 
phenomenology of the high T c systems. 

II. CONSTRUCTION OF PROJECTED SO (5) MODELS 

We begin with the symmetric 50(5) Hamiltonian defined on a lattice, 
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H = AY,L 2 ab (x)- J n a (x)n a {x') 

x <xx'> 

+ V L ab {x)L ab (x') 



(2.1) 



where n a (x) denotes the five component superspin vector on a given site, and L ab (x) is the 5*0(5) symmetry generator, 

Lab = n a p b - n b p a (2.2) 

expressed here in terms of the superspin vector n a and its canonically conjugate momenta p a , 

[n a ,Pb] = ibab (2.3) 

This lattice quantum non-linear a model can be rigorously derived as the low energy limit of a microscopic 5*0(5) 
ladder model [|]|6|. On the ladder, the rung 50(5) singlet state is the vacuum |fi), from which the lowest S0(5) 
£1) is created by a quintet of Bose creation operators, which satisfy 



multiplet | a) = t J 



[t a ,tl] = $ab, t a \n) = o 



(2.4) 



Here, a = 2, 3, 4 denote the triplet (magnon) states, and a = 1, 5 are the hole and particle pair states, (see Fig. ([!])) 
The superspin coordinates are microscopically constructed using these lattice bosons, 



(2.5) 



ta\Q> 
t n |H> 
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FIG. 1. Schematic representation of the singlet state, the triplet magnon states and the hole and particle pair states. 



Due to their microscopic origin, these bosonic states a re h ard-core bosons, in the sense that one can not define 
two of them on the same rung. The A term in equation (|2.l[) describes the gap energy of the magnon and the pair 
states, the J term stands for the hopping and the spontaneous creation/destruction process of these states, and the V 
term describes their nearest-neighbor interaction. This quantum nonlinear a model can in principle also be derived in 
higher dimensions from a microscopic 50(5) symmetric model by introducing a superspin vector as a Hubbard- 

Stratonovich decoupling field, and integrate out the fermionic degrees of freedom in a gradient expansion. However, 
we shall proceed more heuristically here [p5| . For a two dimensional system, one can imagine that the quantum a 
model Hamiltonian is obtained from a "block spin" type of coarse-graining of the microscopic electron Hamiltonian, 
and is defined on a lattice with twice the lattice constant compared to the microscopic electron model. (This doubled 
unit cell is the minimal size needed to define the local AF and d wave SC order parameters). Therefore, each site x 
in the effective model correspond to a plaquette of the microscopic electron model. (On a ladder, this corresponds 
to going from the lattice sites to ladder rungs). The 50(5) singlet state |f2) corresponds to a "RVB" type of singlet 
state, while the five-fold states t a \£l) describe the triplet magnon states, and the d-wave hole and particle pair states 
on a plaquette. Unlike the ladder case, the magnon and the d- wave pair states could condense in the ground state to 
form AF and SC broken symmetry states. In fact, Eder [^] has recently shown that properties of the AF states can 
be described by a coherent state of magnon condensation on top of a uniform spin liquid state. Our model therefore 
describes competition among the "RVB" type of singlet vacuum and the two forms of broken symmetry order. 

While it is reasonable to take J and V to be approximately equal for magnons and pairs, the gap energy A for the 
neutral magnons and the charged pairs are very different in the insulating state at half-filling. In fact, their difference 
is of the order of the insulating gap U at half-filling. Taking into account the hard-core condition and neglecting the 
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nearest- neighbor interaction V for now (it has higher powers of time and space derivatives in the continuum limit), 
we can express the general anisotropic 50(5) model as 

H = A s tlt a (x) +A C ^2 Au{x) 

X X 

- J s 2_j n a(x)n a (x') - J c 2J rii{x)m{x') (2.6) 

<xx'> <xx'> 

In this paper we shall use the convention where a,b,.. — 1,2,3,4,5 denote the superspin indices, a,/3,.. — 2,3,4 
denote the spin indices and i,j = 1,5 denote the charge indices, and repeated indices are summed over. The main 
focus of our paper is to consider the limit where A c > > A s . 
Let us define the charge eigen-operators and t p as 

h = ^(t h + t p ) t B = jj=(t h -t p ) (2.7) 

From this definition, it is clear that t' h is the creation operator for a hole pair and t p is the creation operator for a 
particle pair. We can introduce a chemical potential term 

H»=^(4t p (x)-tlt h (x)) (2.8) 

to describe the effects of doping. In the presence of this chemical potential term, the gap energy of the hole and 
particle pairs are A c — [i and A c + fj, respectively. A chemical potential of the order of the charge gap A c is needed 
to induce a metal-insulator transition in this system. Near such a transition point, the gap energy of the hole pair 

A C = A C - M (2.9) 

can be comparable to the spin gap A s , while the gap towards a particle pair excitation is of the order of twice the 
charge gap, and needs to be projected out of the spectrum in the low energy limit. 

Therefore, within this formalism, the Gutzwiller projection is equivalent to restricting ourselves to the projected 
Hilbert space where 

«p(a?)]*> = (2.10) 
at every site x. Within this projected Hilbert space, the projected Hamiltonian takes the form 

H = A s t^t a (x) + A c m {x)rii (x) 

X X 

- J s ^2 n a (x)n a (x') — J c ^2 ni{x)ni{x') (2-11) 

<n'> <n'> 

This Hamiltonian has no parameters of the order of U, and it is reasonable to expect A s ~ A c and J s ~ J c . We 
see that the form of the Hamiltonian hardly changes from the unprojected model, but the definition of n\ and ti 5 is 
changed from 

71! = ^={h + t\) = ^{t h + t p + 4 + *t) 

n 5 = ^h + 4) = - ^ - 4 + 4) ( 2 - 12 ) 



to 

m 



\{th + t{) n 5 = ^(i h -4) (2.13) 



From equation (2.12), we see that ri\ and 125 commute with each other before the projection. However, after the 
projection, they acquire a nontrivial commutation relation, as can be seen from equation ( [2.13| ): 

[m,n 5 }=i/2 (2.14) 

Therefore, the Gutzwiller projection can be analytically implemented in the SO(5) theory by retaining the form of 
the Hamiltonian and change only the quantization condition. 



4 



III. ANALOGY WITH LOWEST-LANDAU-LEVEL PROJECTION 



The discussions outlined above reveal a deep connection between the Gutzwiller projection within the 50(5) 
formalism and the projection onto the lowest Landau level (LLL) in the context of the fractional quantum Hall 
effect. Consider the problem of a charged particle in a strong magnetic field B and a rotationally symmetric potential 
V(X,Y). In the absence of a magnetic field, all eigenstates form irreducible representations of the two dimensional 
rotation group 0(2), characterized by integral eigenvalues of the angular momentum operator 

L Z = XP Y -YP X (3.1) 

However, in the presence of a strong magnetic field and projected into the LLL, only negative eigenvalues of Lz are 
realized. This is analogous to the situation encountered here. The local charge operator in the 50(5) theory takes 
the form of the angular momentum in the ri\ — n§ plane as given by equation ([0]). When doubly occupied sites are 
locally projected out, the local charge operator, or the angular momentum in the n\ — n§ plane, takes only negative 
values. Since the chemical potential couples directly the angular momentum in the n\ — n<^ plane, it plays the role 
of a fictitious magnetic field threading every £0(5) rotor in the ri\ — 715 plane. The Landau level spacing Tiuj c is 
analogous to the charge gap A c encountered here, and both are taken to be infinity in the projected models. After 
the projection, the Hamiltonian in the Landau level problem retains its 0(2) symmetric form, 

H = V(X,Y) (3.2) 



although a new quantization condition is imposed between X and Y, as given by equation (1.3). This is analogous 
to the observation we made here that the Hamiltonian formally retains a SO (5) symmetric form after the projection 
( [2.10 ), but the quantum dynamics is changed due to the non- trivial commutator between n\ and n^. In both cases only 



a part of the full symmetry multiplets remain after the projection. However, the formal symmetry of the Hamiltonian 
has direct physical manifestations despite the projection. For example, in the LLL problems, semi-classical orbits of 
the guiding center coordinates are still 0(2) symmetric. In our case, we shall see that the static potential for the 
superspin vector can still be 50(5) invariant despite the projection. 

Perhaps the most explicit way to establish the precise connections between these two problems is to consider the 
constraints on the wave function. In the symmetric gauge of the LLL problem, the annihilation operator for the 
cyclotron coordinates takes the form 

a^dg + z/A (3.3) 
where z = X + iY and z = X — iY. Projection onto LLL requires 

a%(z,z) = (3.4) 

which determines the form of the LLL wave function to be 

*(z,z) = /(z)e~^ 4 (3.5) 

where f(z) is a holomorphic function of z only. This holomorphic condition also places strong constraints in many- 
body systems and led to the celebrated Laughlin's wave function. Our no-double-occupancy constraint ( 2.1Cf ) is 



analogous to the LLL constraint (3.4). In fact from equations (2.7) and (2.12), we obtain 



t P = \{z + 2dz) (3.6) 

where z = n± — in§ and z = n\ + in§. For a single unprojected 50(5 ) rotor, the wave function ^(n a ) is a function of 
the superspin coordinates. However, the Gutzwiller projection ( |2.10| ) restricts the wave function to be 

*(ni, n 2 , n 3 , n 4 , n 5 ) = f(z = n\ - in 5 ,n 2 ,n 3 ,n i )e^ zz/2 (3.7) 

where f(z, n 2 , n-3, 714) is a holomorphic function of z. For a collection of 50(5) rotors, the superspin coordinates are 
themselves functions of the lattice sites x, and \P[n a (x)] is a functional is the superspin coordinates at each site. For 
the projected 50(5) models, this functional is restricted to take the form 

*K(.x)] = f(z(x), n a {x)) J] e -^(-)/ 2 (3.8) 
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where f(z(x), n a (x)) is a holomorphic functional of z(x) — ni(x) — in$(x). 

The formal but precise analogy between the two types of projection allows us to introduce the concept of a chiral 
50(5) rotor. This is a system of rotors with 50(5) invariant potential and coupling, however, the rotation within the 
ni — 715 plane is chiral, i.e. only one sense of the rotation is allowed. Such a system of chiral 50(5) rotors is described 
by the wave functional in equation (|3.8|). 



IV. SO (5) SYMMETRY OF THE GROUND STATE ENERGY 

Having discussed the general notions of the projected 50(5) model, we are now in a position to explore the phase 
diagram of this model. As we commented earlier, the projected 5*0(5) model describes the competition and unification 
of the spin liquid, AF and the SC states. In the original unprojected 5*0(5) symmetric model, not only are the AF and 
SC states degenerate in energy, but they are also degenerate with all the intermediate coexistence states. This points 
out a route from AF to SC with no potential barrier, and introduces the concept that the metal-insulator transition 
in the high T c systems can be viewed as a smooth rotation of the 50(5) superspin. One of the key questions to be 
answered in this work is what happens to the picture in the case of projected 50(5) symmetry. 

In anticipation of the competition of the states discussed above, we construct a class of variational wave functions 
in the coherent state representation: 

l*> = Y[(cos6(x) +sm9(x)(m a (x)ti(x) + A(x)t{(x)))\n) (4.1) 

Here |fi(x)) denotes a local singlet state defined by t a (x)\£l(x)) — th(x)\fl(x)) = and |0) is a product state of these 
local singlets, |fi) = J\ \&( x ))- S(x) is a local variational parameter describing the competition between long range 
order and quantum disorder. For 0(x) = our variational wave function describe a spin singlet ground state, while 
a non-zero value of 9 describes a coherent state formed by the local singlet, the magnon or the hole pair state. This 
wave function is a generalization of the coherent state description of a AF state in terms of a magnon condensate 



5j. As we shall see from the following equation (4.2), sin29 stands for the length of the 50(5) superspin vector. 



m a (x) and A(x) are general complex variational parameters describing the local amplitude for magnons and hole 



pairs. We notice that this wave function satisfies both the Gutzwiller constraint (2.10) and the hard-core constraint 
for magnons and hole pairs exactly. It is easy to see that 

(*|n a (a;)|*) = -^=sin29(x)Re(m a (x)) 
v2 

(#|ni(a;)|#) = -sin29(x)Re(A(x)) 

(*|n 5 (x)|*) = ~sin29(x)Im(A(x)) (4.2) 

where Re and Im denote the real and imaginary parts of a complex number. The coupling terms in the projected 
50(5) Hamiltonian depend only on n a (x), ri\(x) and 71.5(2;). Therefore, the coupling energy depends only on the real 
part of m a (x) while it depends on both the real and imaginary parts of A.(x). Therefore, for discussing the ground 
state wave functions, we can assume without loss of generality that m a (x) is real and A(x) = mi(x) + im^(x). The 
normalization condition (^l^) = 1 can be implemented by the constraint that 

m 2 a (x) = (m\ + m\ + m 2 a ){x) = 1 (4.3) 

Therefore, we see that although we have completely projected out the particle pair states, the local degrees of freedom 
can still be represented by a vector on a five dimensional sphere. 

Uniform states are obtained by taking all parameters to be constant. For A = and sin29 ^ 0, our wave function 
\ty) describes a pure AF state with the following properties: 



(¥|Q|¥> = <*| 5>J,*fc(:0l*> = 

X 

(*|JV a |*) = (*!]>>„ (X) I*) =N^-=sin29 m a 

X * 
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*|5 2 |*) = (*|(^5(a;)) 2 |*) ^N2sm 2 9 + Nsm i e(l~m 4 a ) (4.4) 



where N is the number of lattice sites. Equation (4.4) describes a half-filled state with a macroscopic Neel mag- 
netization, and vanishing uniform magnetization. Furthermore, this state is composed as a linear superposition of 
eigenstates with different values of the total spin, and the fluctuation of the total spin scales like J (S 2 ) oc y/N, just 
as one expects from a standard Neel state. 

On the other hand, for m a — and sin20 ^ 0, |\&) describes a pure SC state with the following properties: 

(*|Q|*> = (^\J2 t h t h(x)m = N S m 2 e 

(*|JVi + iN 5 \\V) = (*| ^2(ni + in 5 )(x)\\V) = N^sin29 (m x + im 5 ) 

X 

(*|Q 2 |*) - (I'M*) 2 = N sin 2 9 cos 2 9 (4.5) 

Equation ( (4.5| ) describes a state with a finite doping density and a finite SC order parameter. Just as in the standard 
BCS case, this state is composed as a linear superposition of eigenstates with different values of the total charge, and 
the fluctuation of the total charge scales like \J (Q 2 ) — (Q) 2 oc V^V, just as one expects from a standard SC state. 

However, besides these two pure states, there is a class of mixed states which interpolates between the pure AF and 
SC states. Taking mi = sin a and = cos a, we see that the mixed states have the following property: 

(#|Q|#) = N sin 2 9 sin 2 a 

(*|JVi +iN 5 \V) = N-sin26sma 

(*|JV 2 |*) = N—=sin20cosa (4.6) 
v2 

Therefore, we see that there is a continuous family of intermediate mixed states interpolating between the pure AF 
state at half-filling and the pure SC state with finite doping density. As the SO (5) angle a rotates continuously 
from a pure AF state with a = to a pure SC state with a — tt/2, the hole density of the mixed state interpolates 
continuously between these two limits. Therefore, our wave function gives a unified description of AF and SC and 
points out a precise route from AF to SC as the doping level is varied. In order for this route, or small deviations 
from this route, to be physically realized in the high T c superconductors, we have to demonstrate that there is no 
large energy barrier for the intermediate mixed states, or that the ground state energy is approximately independent 
of the 50(5) mixing angle a. In particular, we have to show that the energy barrier is independent of the Hubbard 
energy U, in the limit of large U. In the following, we shall investigate this question. 

The energy functional |W) describes the coupling between these five dimensional vectors ra (i), and it is given 

by 

(\V\H\V) = E(9(x),m a (x)) 
= — — sin2O{x)sin20(x f )in a (x)i7i a (x f ') — sin26(x)sin29{x f )mi(x)rrii{x r ) 

xx' xx' 

+A S ^ sin 2 6»(x)m 2 (a;) + h c ^s\n 2 9(x)m 2 (x) (4.7) 



This ground state energy functional describes a systems of coupled rotors satisfying the constraint (4.3). At the point 
J c = 2J S and A c = A s in parameter space, this rotor model is exactly 50(5) symmetric. This is a central observation 
of this work. 

From this consideration we learn a very important lesson about the compatibility of the Mott insulating gap and 
the idea of a smooth 50(5) rotation from AF to SC. As we have seen, the large asymmetry between the charge and 
spin gap at half-filling necessitates the removal of the Q(x) > part of the 50(5) multiplets, therefore, the dynamics 
close to half-filling has to be modified. But the static potential governing the transition from AF to SC can remain 
50(5) symmetric, and in particular, the energy barrier separating these two states can remain small in the limit 
where the Mott insulating gap tends to infinity. 

We also observe a crucial difference between the projected and the unprojected 50(5) models. In the unprojected 
50(5) model with the full 50(5) symmetry, AF and SC states are degenerate at half-filling, and the rotation between 
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these two states can be continuously performed without changing the density to going away from half-filling. This 
case is similar to the well-known degeneracy between the CDW state and the s wave SC state for the negative U 
Hubbard model at half-filling. In the projected 50(5) model, where all particle pair states have been locally removed, 
the 50(5) rotation from the AF to SC states are accompanied by the continuous change of the hole density, and a 
pure SC state can only be reached at a finite critical hole doping density p c — sin 2 8. While the unprojected 50(5) 
symmetry is only valid at half-filling, the projected 50(5) symmetry can be valid for a range of doping concentration 
< p < p c , since all these doping concentrations correspond to the same value of the chemical potential p = p c at 
which the ground state energy (4.7) is S0(5) symmetric. 

However, it should be pointed out that the projected 50(5) symmetry at p = p c has only been demonstrated 
within the variational mean field approximation. This corresponds to the semiclassical limit, and becomes exact only 
in the large s limit, where s labels the representation of the local 50(5) group at a given site. Quantum fluctuations 
can be systematically investigate d as a 1/s expansion. Assuming uniform ground states, we have studied the effect 
of zero point fluctuations in Sec. VII and found that at p, = p c , the intermediate mixed states have slightly higher 
energy than the AF and SC state. Therefore, quantum fluctuation leads to a slight breaking of the projected 50(5) 
symmetry. The important point here is that this symmetry breaking effect can be systematically controlled in the 
semi-classical 1/s expansion, and certainly is independent of the Coulomb energy scale U. This fluctuation would 
induce a first order transition and predict phase separation of AF and SC states at p = p c . However, there are also 
other competing interactions such as nearest-neighbor and next-nearest-neighbor interactions which tend to reduce 
the barrier, and could also lead to non-uniform states like stripes. Due to the complexity of the calculations, we shall 
defer the detailed studies of these competing effects to future works. 



V. PHASE DIAGRAM 

In this section we investigate the phase diagram of the projected 50(5) model within the framework of the variational 



wave function (4.1). Taking uniform values of the variational parameters 8, m x = cos a, mi = sin a and m y 



m. 



TO5 = 0, the variational energy in (4.7) reduces to 



E (^ a ) = _ j s s i n 2 28 cos 2 a- y sin 2 28 sin 2 a + A s sin 2 8 cos 2 a + A c sin 2 8 sin 2 a (5.1) 

,2 , 



In the following, we shall mainly study the 50(5) symmetric case, and take J c = 2J S = 2J. Defining x — sin 8 and 

Efeg) x — A& ar ,H X = As. 

iJN ' ° s - 4J ana ° c - 4J • 

e(x, y) = x 2 - x + 5 c x + (5 S - 5 c )xy (5.2) 



y = cos 2 a, and the dimensionless coupling constants e = E }^'^ , 5 S = xf an d § c = Tf> we obtain 



We shall minimize (5.2) with respect to x and y, subject to the condition that < x, y < 1. The phase diagram 
can be plotted in the two dimensional parameter space of S c and S s . We notice that e(x,y) depends linearly on y, 
therefore, for S s > 5 C we obtain y m in — and 

< x min = < 1 for -1<S C <1 

Xmm = for 8 C > 1 

X m in = 1 for S c < -1 (5.3) 
Similarly, for 6 S < S c we obtain y m i n = 1 and 

< x mm = < 1 for -1<8 S <1 

Xmin = for S s > 1 

Xmin = 1 for S s < -1 (5.4) 

From these equation, we can determine the phase diagram as shown in Fig. (0). 
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0<x<l, y=0 
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x=l, y=l 
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FIG. 2. Phase diagram of the projected SO(5) model in the S s versus 8 C plane. Phase boundaries are depicted by the solid 
lines. Variation of the chemical potential traces out a one dimensional trajectory as shown on the dotted line. 

There are seven different phases on this phase diagram. x m i n = corresponds to a quantum disordered singlet 
state with no condensed bosons. x m in = 1 and y m in = corresponds to a quantum disordered state with completely 
filled hole pairs. x m in = 1 and y m in — 1 corresponds to a quantum disordered state with completely filled magnons. 
< x m i n < 1 and y m in = 1 describe a pure AF phase, while < x m in < 1 and y m in = describes a pure SC phase. 
When — 1 < 6 S = S c < 1 a continuous family of mixed AF/SC states labeled by a free superspin angle < a < ir/2 is 
realized, while for S s = 5 C < — 1 a continuous family of quantum disordered states labeled by a free superspin angle 
0<a<7r/2is obtained. 

The system traces out a one dimensional trajectory in this two dimensional phase diagram as the chemical potential 
is increased, as depicted in Fig. (||). Increasing the chemical potential decreases the S c parameter while holding 8 S 
constant. S s describes the degree of quantum spin fluctuations in the system, since in the AF phase, the size of the 
Neel moment 



decreases with increasing 5 S . For S s < 1, the system goes through a phase transition from AF to SC at = /i c = 
A c — A s . At this critical value of the chemical potential, the 9 parameter remains fixed, but the a parameter changes 
continuously from to 7r/2, and correspondingly, the density changes from to 




(5.5) 



_ - 2 n _ _ L "s 

Pc — sin (7 x m i n ^ 
This behavior gives a density versus chemical potential diagram as shown in Fig. 



(5.6) 
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He h 

FIG. 3. Density versus chemical potential relation in the projected SO(5) model. Unlike the case of a generic first order 
transition, the ground state in the density range < p < p c is uniform, rather than phase separated. 

As we see, for densities in the range < p < p Cl the system is infinitely compressible since dp/ dp, = oo. For // > p c 
or p > p Cl the system has a finite compressibility of dp/dp = 1/8J. 

It is interesting to plot both the AF and the SC order parameters (ni) and (712) as a function of the density for 
the whole range of < p < 1. We will restrict to the case of < S s < 1 where the undoped state is a AF state. We 
obtain the following doping dependence of the SC order parameter: 



(m) = $ V 7 /^ 1 - Pc) for p< p c 



vW - P) f or P> Pc 
and the doping dependence of the of the AF order parameter: 

/„ \ _ / V 2 ( l - Pc){Pc - P) for p< p c 
{n2) ~\0 for p>p c 

These behaviors are depicted in Fig. (|J). 



(5.7) 



(5.8) 



M Order Parameter 




\ <n 2 > 






\ <ni> 



p c 0.5 1 P 

FIG. 4. AF ((712)) and SC ((wi)) order parameters versus density in the projected SO(5) model. 



We note several interesting features of the phase diagram. First of all, we can use the energy functional (4.7) as 
a starting point for a finite temperature classical fluctuation analysis, and estimate the transition temperature due 
to the classical fluctuations. Within such a framework, the 3D AF (Tjv) and SC (T c ) transition temperatures are 
proportional to the stiffness of the spin and the phase fluctuations, which are in turn proportional to (ri2} 2 and (rii) 2 
respectively. Therefore, Fig. (^) gives an approximate estimate of the transition temperatures. For small p c , we see 
that there is a sharp drop of Tjv and a maximum of T c at p = 1/2. The system is a AF insulator at doping p = and 
a pure SC state for p > p c . The reason for a maximum of T c at p — 1/2 is due to the strong correlation of the charged 
bosons. The charge bosons have a hard-core interaction, therefore, they are insulating at both p = and p = 1 and 
have maximal charge stiffness at p — 1/2. We can perform a rough translation of this optimal doping value in our 
effective model to the microscopic model. Since our effective model is defined on a unit cell with twice the lattice 
spacing of the microscopic model, p = 1/2 therefore describes one hole pair per four sites in the microscopic model, 
or a doping of x = 1/4 = 25% in the conventional language. This crude argument tends to overestimate the value for 
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optimal doping, since it neglects the effects of unpaired electrons. However, considering the crudeness of the estimate, 
it is still reasonably close to the optimal doping x = 15% observed in the LSCO family of high T c superconductors. 

In the regime of < p < p c , the system is a coherent mixture of AF and SC order. For this entire range of densities, 
the system has a projected 50(5) symmetry within the variational approximation discussed above. The projected 
50(5) symmetry manifests itself in terms of a infinite compressibility in the region < p < p c and, as we shall see 
in next section, a charge mode with a dispersion relation lu ~ k 2 . Since such a state is rather unusual, and maybe 
highly susceptible to density fluctuations, we would like to discuss more detailed physical properties in this region. 

First let us comment on the fact that there are familiar physical systems whose uniform ground states are infinitely 
compressible. The free boson model is certainly such an example, and the density mode also has a uj ~ k 2 dispersion 
relation. But the infinite compressibility is due to the absence of the interaction, which is not characteristic of the 
strongly interaction system considered here. A less trivial example is the spin 1/2 XXZ ferromagnetic Heisenberg 
model, given by the Hamiltonian, 



where J < and the sum extends over nearest neighbor sites of a square lattice. This model can be interpreted as 
quantum hard-core boson model, where the fully polarized spin down state could be identified with the vacuum of the 
bosons, the XY part of the Hamiltonian describes the hopping of the bosons and the last term describes the nearest 
neighbor attraction between the bosons if A > 0. When A > 1, the system is in the Ising limit, and the spontaneous 
breaking of the Z2 symmetry implies phase separation of the bosons. On the other hand, when < A < 1, the system 
is in the XY limit, and the ground state is a superfluid. Therefore, the anisotropy parameter describe the competition 
between superfluidity and phase separation. At A = 1, the system has a SU(2) symmetry and the dispersion relation 
becomes quadratic. Different directions of the ferromagnetic polarizations are degenerate and can be changed without 
any energy cost. Since the z component of the ferromagnetic polarization is identified with the total density of the 
hard-core bosons, the SU (2) vacuum degeneracy implies infinite compressibility of the corresponding boson system 
for the entire range of boson densities < p < 1. 

These two examples illustrates that there is nothing intrinsically pathological about having a system with infinite 
compressibility. The second example is more generic, and shows that uniform states with infinite compressibility 
can be obtained in systems on the verge of phase separation, and the infinite compressibility can be ensured by 
symmetry. Both of these properties are also shared by the projected 50(5) model. These models are on the verge 
of phase separation into AF and SC phases, and the infinite compressibility is a result of the projected 50(5) 
symmetry. It indicates that small perturbations, such as next nearest neighbor interactions, quantum fluctuations, 
and quenched disorder will be very important to determine the true ground state. With such perturbations, the 
ground state is expected to be unstable toward the experimentally reported textures such as the spin glass, stripes 
and incommensurate spin density waves. 

Next let us investigate the phcnomcnological consequence of this remarkable property. One of the most puzzling 
properties of the high T c superconductors is the constant chemical potential in the underdoped samples. For LSCO 
systems, where the doping level can be varied continuously by the 5r concentration, this effect has been dramatically 
observed in the ARPES experiments and the constant chemical potential persists from the weakly doped insulator 
to optimally doped superconductor pq ). In fact numerical calculations on the Hubbard model also reveal similar 
divergent behavior of the compressibility as the metal- insulator transition is approached from the metallic side p7[ | . 

The simplest explanation of the small chemical potential shift is a two phase mixture with different densities at 
a first order phase transition. If the system globally phase separates into two different spatial regions with different 
charge densities but the same free energy densities, the added charges only change the proportion of mixture of the 
two phases and do not change the energy, therefore, dp/ dp = 0. However, this situation of global phase separation 
can certainly not occur in a system with long ranged Coulomb interaction and is ruled out in the real high T c system. 

A phenomenon possibly related to the tendency of phase separation is the formation of stripes [2g|-|3C|l . A stripe state 
can be viewed as microscopic phase separation of AF and SC into alternating regions, where each region has different 
charge density and the same free energy density. However, a crucial difference between the global phase separation 
and this picture of microscopic phase separation is that the stripe state has infinitely many surfaces between AF and 
SC, and the surface energy makes a finite contribution to the total energy in the thermodynamic limit. In this picture, 
doping can be accomplished by converting AF stripes into SC stripes, thereby creating more surfaces separating AF 
and SC regions. Finite doping density therefore leads to a finite density of surfaces and the accumulated surface energy 
would in general lead to a shift of the chemical potential under a generic situation. Additional physical conditions 
are needed to ensure the constant chemical potential in the stripe phase. 




(5.9) 
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Therefore, the absence of the chemical potential shift places a very strong constraint on possible theoretical expla- 
nations. The projected 50(5) model proposed in this work offers a possible explanation for the absence of chemical 
potential shift. At a critical value of the chemical potential where the AF and SC states have degenerate energy 
density, we can have three situations, where the intermediate mixed states have higher, lower or degenerate energy 
compared to the AF and SC states. When the intermediate states have higher energy, the system will go through a 
hrst order phase transition at p, = p, c and this will lead to global phase separation into AF and SC regions. On the 
other hand, if the intermediate states have lower energy, there exists a range of chemical potential p c \ < p, < p c % 
where the mixed phase has a uniform and continuously varying density. In this case, dp/dp ^ is obtained. When 
the region p c \ < p < p c i shrinks to zero, we obtain the limiting 50(5) symmetric case where the system is on the 
boundary between a first order transition and two second order phase transitions. In this case, dp/ dp = for a range 
of densities < p < p c . Therefore, if we restrict to ground states where the density is not globally inhomogeneous, 
the absence of the chemical shift directly implies 50(5) symmetry. 

Within this model, we can therefore define a experimental procedure to measure one of the most crucial parameter of 
the theory, namely p c . Using the experimental ARPES data for LSCO system we would identify p c to be approximately 
the same as the optimal doping density. We recall that p c is also a measure of the degree of the quantum spin fluctuation 
in the system. For the bi-layer materials such as YBCO and BISCO superconductors, the quantum spin fluctuation 
are stronger due to the inter-layer spin exchange, and we would predict that p c should be less than the optimal doping 
value. 



VI. COLLECTIVE MODES 



Having discussed the ground state properties and the phase diagram of the model, we are now in a position to 
study the collective excitations of the model. We have argued that the ground state energy can remain 50(5) 
symmetric despite the projection. However, the projection does affect the collective excitation spectrum near half- 
filling. Nonetheless, as we shall see, there remains a unique signature of the projected 50(5) symmetry in the collective 
excitation spectra. 

In principle, the collective excitation spectra can be obtained straightforwardly by studying the quadratic fluctu- 
ations around the mean field minima. The resulting quadratic boson Hamiltonian can be simply diagonalized. The 
main complication in the procedure is the hard-core boson constraint, which requires 

tit a (x)+t{t h (x)<l (6.1) 

for every site. There are several ways to implement this constraint rigorously. One is to follow the mapping from 
the one-component hard-core boson model to the XY model and generalize it to a multi-component hard-core boson 
model. One could also convert the above inequality constraint to an equality constraint by introducing a boson 
creation and annihilation operator for the singlet state. This approach will be implemented in the Appendix. For 
simplicity of presentation, here we shall adopt a less rigorous approach and introduce a on-site boson repulsion term 

Wj2(tU<*+ttt h ) 2 (6.2) 



to our Hamiltonian ( |2.11 ) and convert the hard-core constraint to a soft-core constraint. We shall show later that 



all results can be expressed in terms of the order parameter, which is implicitly dependent on W, but there is no 
explicit dependence on W. We have verified that all three methods give the same long wave length spectra for the 
collective modes in the limit of low boson density. In the next section, we shall present another calculation based on 
the continuum effective Lagrangian method, which also reproduces the same spectra. 



To simplify presentation, we shall concentrate on the case where the ground state energy functional (4.7) is 50(5) 
symmetric, i.e. for coupling constants J c = 2J S = 2 J and A c = A s = A. We choose the direction of spontaneously 
broken symmetry to be (t x ) = (t x ) = x and (th) = (ty = y. The extremal condition can be easily determined to be 

^ +yW = f^ (6 . 3) 

The combination x 2 + y 2 expresses the fact that the classical minimum is 50(5) symmetric. We can therefore write 
x = r cos a and y = rsina. Expanding the boson operators as: 

t x =x + a x , t y = a v , t z = a z , t h = y + a h (6.4) 
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we obtain the following quadratic Hamiltonian 



H = (A + 2Wr 2 ) ^(4 a x + a{a h ) + Wx 2 ^(4 + a,) 2 + Wy 2 ^(a* + a hf 

X XX 

- 2 + 51 (flh{x)a h {x') + h.c.) 

<x,x'> <x,x'> 

+ 2WxyY,(a{(x) + a h (x))(al(x) + a x (x)) 

X 

~2 X! K( a; ) + a a( a; ))K( a; ') + a a (^))- g I] (4(z) + a*(*))(4(z') + a*(z')) 

<a: ) .a: / > <sc ) a: / > 

+ (A + 2Tyr 2 )^(4 a , + 4 az ) (6.5) 

X 

We are in particular interested in the collective mode spectra for the AF insulating state with a — 0, the mixed states 
with < a < tt/2 and the SC state with a = tt/2 and how they connect to each other. 

From this quadratic Hamiltonian we can learn a number of important features. First we notice that the a y and a z 
modes are decoupled for all ranges of < a < tt/2, but most importantly, their dispersion relations are independent 
of a and given by 

uj(k) = v s k , v s = 2 J (6-6) 

where k = a\k\ and a is the lattice constant. This is indeed a very remarkable property. At a — 0, a y and a z modes 
are nothing but the transverse AF spin wave modes. AF spin waves are usually viewed as Goldstone modes and their 
existence is due to the AF long range order. However, as a changes continuously from to tt/2, the AF long range 
order continuously diminishes until it vanishes at a = tt/2. The reason that the properties of the a y and a z modes do 
not change at all is due to the SO(5) symmetry of this model, since the diminishing AF order is compensated by the 
increasing SC order as the superspin angle a is varied. As we shall see, the a x mode is the AF spin amplitude mode 
at a = 0, but it becomes massless and degenerate with the a y and a z modes at a = tt/2. These three modes form 
a massless tt triplet mode whose existence is purely a consequence of the SC order. Therefore, as a is continuously 
varied from to tt/2, the transverse AF spin wave modes gradually change their character to become the tt triplet 
resonance of the SC state. As we shall see, for /i > fi c , the tt triplet mode becomes massive. 

At the AF point a = 0, the spin amplitude mode a x is decoupled from the SC mode ah and can be diagonalized 
separately. The dispersion for the spin amplitude mode has the conventional massive relativistic form for small k 

^(k) = 16J 2 (fc 2 /4 + m 2 x ) , m 2 x = ™f = (6.7) 

On the other hand, we have a massless SC Goldstone mode ah with the following dispersion, 

u){k) = Jk 2 (6.8) 

This mode is an important new prediction of the SO (5) theory. It is the counterpart of the tt resonance in the 
AF state. In the unprojected 50(5) model, there are two such modes, with charge ±2, and they represent gapless 
fluctuations from AF to SC at half-filling. In the projected SO (5) model, the charge +2 mode is projected out of 
the spectrum, however, the charge —2 mode remain massless at = /i c . It is also a manifestation of the gapless 
fluctuation from AF to SC at half-filling, but the SC fluctuation is hole-like, rather than both hole and particle-like as 
in the unprojected case. We see again that a large Mott-Hubbard gap is fully compatible with gapless SC fluctuation 
at half-filling. Experimental detection of this mode could provide a important test of the projected SO(b) symmetry. 
The a x and the ah modes also decouple in the pure SC state with a = tt/2. However, their physical interpretation 



change. The a x mode becomes gapless at this point with the same dispersion as in (3.6). Therefore, the three modes 



a x , a y and a z form a gapless tt triplet mode of the pure SC state, and represent the gapless fluctuation from SC to 



AF at fi — fi c , but with a finite hole density p c given in (5.6). The dispersion for the ah mode is given by 



co(k) = 2Jm h k , m l = ^l = ^L-A (6 .9) 
and has the natural interpretation of a linearly dispersing phase Goldstone mode of the SC state. 
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In the intermediate mixed state with < a < it/2, the a x and the the ah modes are coupled. Diagonalization of 
these modes gives: 



QJ 2 {k) _ \ (l + m 2 )fc 2 /4 + m 2 f m 2 = cos 2 a 

(4J) 2 ~ \ (1 + m 2 h /m 2 x )k A /16 { m 2 = ^-^ sin 2 a 



h> L x 1 ""x — /2j . „" (6 10) 

27 



The upper massive mode has predominantly spin amplitude character, and we see that the gap diminishes continuously 
until it reaches zero at a — tt/2 to become the massless tt triplet mode. The lower mode has predominantly SC 
fluctuation character, and has a gapless u> oc k 2 dispersion. In the mixed region where both the AF and the SC order 
parameters are non-zero, one would naturally expect a gapless phase mode corresponding to the SC order. However, 
in a interacting boson system, the phase mode is expected to have linear dispersion on general ground. Therefore, 
what is surpassing and new here is not the gapless nature of the SC mode, but its quadratic dispersion. In order to 
locate the origin of the quadratic dispersion, we have perturbed the model away from the projected 5*0(5) symmetric 
point so that a uniform mixed state is stabilized as a classical minimum. SC fluctuation around such a non-50(5) 
symmetric point is gapless and has linear dispersion. A quadratic dispersion is only realized at the 5*0(5) symmetric 
point. Therefore, the quadratic dispersion is a unique signature of the projected 50(5) symmetry in the entire range 
of densities < p < p c \ To understand the physical origin of this remarkable phenomenon, we notice that a boson 
system with gapless quadratic dispersion generally has infinite compressibility. This can be directly seen from the 
compressibility sum rule 

K = ^ = i lim lim x(fc,w) ~ lim \- (6-11) 

where x(k,ui) is the dynamical density correlation function. Because of the quadratic dispersion relation, we can see 
explicitly that x(fc) oc 1/fc 2 for small fc, therefore, a infinite compressibility is obtained. On the other hand, a infinite 
compressibility implies that ^ = 0, i.e. the chemical potential is independent of doping. But this is exactly the 
prediction of the projected 50(5) model! For < p < p c , the chemical potential is pinned at the 50(5) symmetric 
point p — pc, where the superspin vector can point in any direction. To accommodate a long wave length fluctuation 
of the hole density, the system rotates into another degenerate minimum with a different superspin angle a and a 
different hole density. For this reason, the chemical potential does not change and the system is infinitely compressible. 



VII. QUANTUM CORRECTIONS TO THE MEAN-FIELD SOLUTION 



At the mean-field level, the ground-state energy of the Hamiltonian ( 2.11 ) with (6.2) depends on the AF and SC 
order parameters x and y only via their combination r 2 = x 2 + y 2 reflecting the SO (5) invariance of the mean- field 
result. However, the zero-point energy of the bosons in the quadratic Hamiltonian gives a correction to the ground- 
state energy due to quantum fluctuations. (Calculations of the quantum fluctuation effects at the 50(3) spin-flop 
transition have been studied in ref. pi| . As we will show below, this correction turns out to depend on x 2 and y 2 
separately. 

For the soft-constraint case, E\ can be expressed as 



(7.1) 



where uj i (k) are the four collective modes described in Sec. VI (although extended to all values of k in the BZ). In 
a systematic 1/s expansion, this correction to the ground state energy scales like 1/s, and is therefore small in the 
semi-classical limit. Alternatively, E\ can be seen as a small correction to the mean-field energy Eq for small values 
of the parameter e = h^j- = 1 — A 
produces a small SO(5)-symmetry breaking, 

2 sin 2 



■^j. However, contrary to Eq, E\ also depends on the superspin angle a and thus 



(7.1) in the small-e limit. We thus parametrize x = 1 



We now evaluate the a-dependent part of Eq. 
(— 1 < X < 1), differentiate the expression (7.1) with respect to Xj an d expand it up to second order in e. By further 



transforming (cos k x + cos k y )/2 = 1 — q with < q < 2, the derivative of E\ can be expressed as the sum of two terms 



— Ei - E' 1A + E[ B , 



(7.2) 



where 
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] A = 2 J f d q V(q) 

IB JO 



G A 
Gb 



with 



° A = ~ 2^— +0{£) 



(7.3) 



(7.4) 



Gb ^ - {V ~ q ~ l] (X/*- 1 -*-^)' 2 + (e3) , (7 .5) 
4 (1 + y/q) 

and the density of states T>(q) is defined as 

T>(q) — — ^ J dk x dky 6 [q — 1 + (cos k x + cos k y )/2] . (7.6) 

Unfortunately, the terms of the e expansion in Gb diverge when integrated over q. It is thus convenient to carry out 
the transformation q = ez in G# and i/ien expand in powers of e. We obtain 

SB = ^f * -±)+0(e). (7.7) 



2 VV^+ 1 +X Vz 



For small e, the integral in z can be extended to oo and we obtain from Eq. (7.3) 

,3/2 



f°° —( /I 4- v f a / 2 l 

^ B = 2 J e / D(«e) G B = VV t X — '- + 0(e 2 log e) . (7.1 

Jo 7T 



While in evaluating Eq. (7.8) we only need the density of states at q — 0, 2?(0) = 7r _1 , for the first term E' 1A one 
needs T>(q) in the whole domain < q < 2 However, is independent of \ an d t nus it merely fixes the value of the 
critical chemical potential. A numerical integration yields 

E' 1A = 0.28 J e . (7.9) 



By integrating over \ Eqs. ( 17-Sh and (7.8), we finally obtain the total contribution to the ground-state energy 
correction. 

J7 = 0.14 X e - A (e (1 + x )f 2 + 0(e 2 loge) + const. . (7.10) 

Ei thus lifts the degeneracy as a function of % and initially favors the pure superconducting phase (x = — 1). A small 

chemical potential term —[i'y 2 = —[i'r 2 (l — x)/2 with ^j- = = —0.28 e + ^ (2e) 3 / 2 restores the degeneracy 
between the pure superconducting (x = —1) and the pure antiferromagnetic (\ = +1) phases. However, due to the 
convexity of E\ as a function of x there is always a barrier (oc Je 3 / 2 ) between the two phases since the mixed phase 
always has a higher energy. This means that at // = (i' c and for intermediate densities the system prefers to phase 
separate between the two pure phases rather than choosing the mixed phase. However, the important point is that 
this barrier remains small in the limit U — > 00. 

In view of the symmetry breaking effe cts of the quantum fluctuations, it would be interesting to see whether there is 



a limit when the the wave function (11) and the projected 50(5) symmetry becomes exact. Rokhsar and Kotliar |32f | 
have shown that this type of wave functions are actually exact in the limit of infinite dimensions. Therefore, besides 
the 1/s expansion, we could also use a l/d expansion (where d is the space dimension) to systematically control the 
50(5) symmetry breaking quantum effects. Besides quantum fluctuations, there are also other symmetry breaking 
terms. Nearest-neighbor and next-nearest-neighbor Coulomb interactions also break the 50(5) symmetry, however, 
their corrections to the ground state is concave, i.e. the energy of the intermediate states are lowered. Therefore, they 
can also lead to uniform mixed states in some region of the phase diagram. The detailed study of all these competing 
effects will be carried out in subsequent works. 
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VIII. LOW ENERGY EFFECTIVE LAGRANGIAN 



While the projected 50(5) model defined on a lattice enables us to make some connection to the underlying 
microscopic physics, for most discussions concerning the long wave length and low energy degrees of freedom, it is 
desirable to have a effective continuum Lagrangian. Such a formulation can be directly obtained by taking the long 
wave length limit of the projected bosonic model discussed previously. However, in order to make the connection to 
the unprojected 50(5) model clearer, we shall motivate our discussion from the original 50(5) effective model. 

The effective Lagrangian for a fully 5*0(5) symmetric model takes the form of 

£ = \{dtn a f - ^(d k n a ) 2 - V{n) (8.1) 

where \ measures the superspin susceptibility, p measures the superspin stiffness, k — x,y denotes the spatial directions 
and V(n) is a scalar function of the superspin magnitude n 2 only. There are three important symmetry breaking 
effects connected with the presence of a large Mott- Hubbard gap. First is an asymmetry in the scalar potential, which 
can be described by a additional term 

V g {n)=+ g -{nl + nl) (8.2) 

which for positive g favors AF at half- filling. Second is the asymmetry between the spin (xs) and the charge (% c ) 
susceptibilities, which modifies the kinetic energy to 

^{d t n a f + ^-{d tni f (8.3) 

The last symmetry breaking effect is due to the chemical potential p, which enters the Lagrangian as a gauge coupling 
in the time direction, and modifies the charge part of the kinetic energy to 

^-{(dt^+^f + idtn^^f) (8.4) 

Combining these three symmetry breaking terms, we obtain 

C = ^(d t n a ) 2 + ^{{d t nx + pn 5 f + (d t n 5 - pn,) 2 ) - ^(d k n a ) 2 - V(n) - |n? 

= ^{dtn a f + ^(d tni ) 2 + fixdnsdtn, - ni d t n 5 ) + ^ Xc ~ 9 n 2 - P -{d k n a f - V{n) (8.5) 

In the presence of a large Mott-Hubbard gap, all these three symmetry breaking terms are of the order of U, i.e. 
Xj 1 ~ 9 ~ A*c ~ U . Therefore, this Lagrangian contains high energy degrees of freedom of the order of U. However, 
as already observed in |lj , at fi = ji c = \J g/xc, their effects cancel completely in the time-independent part of the 
Lagrangian, and the static potential is SO(5) symmetric just as in the original unprojected model. We also observe 
that near the AF/SC transition point where /i ~ fj, c , the first order time derivative term is of the order of one. 
Furthermore, in the spirit of the low frequency and wave vector expansion, we only need to retain the first order time 
derivative term in the charge sector and can drop the second term in the above Lagrangian. Combining all these 
considerations, we obtain the following low energy effective Lagrangian near the AF/SC transition region, which is 
free of any parameters of the order of U: 

C = y (9 t n Q ) 2 + (n 5 d t m - ni d t n 5 ) - ^(d k n a ) 2 - V(n) (8.6) 

This is exactly the Lagrangian counterpart of the £0(5) projection procedure discussed previously in the Hamiltonian 
language. Dropping the second order time derivative terms removes half of the (high energy) degrees of freedom, and 
redefines the canonical conjugacy of the dynamical variables. In particular, the conjugate variable of n\ is nothing but 
715 itself, since p± — 5C/Sdtni = 2n$. Standard quantization procedure requires the canonical commutation relation 
[ni,pi] = i, which in this case just reproduces equation fll.4| ). This confirms the fact that the SO(5) projection does 
not change the form of the interaction potential, only the commutation relation between ri\ and n§. 

It is easy to see that the low energy effective Lagrangian (|8.6|) produces exactly the same long wave length collective 



mode spectrum as the projected S0(5) Hamiltonian (2.11) defined on a lattice. To facilitate the comparison, we take 
the 5*0(5) potential to be 
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(8.7) 



2 ^ " ' 4 

a 

Assuming broken symmetry in the ri\ and n,2 directions, we find that the 77.3 and modes always decouple, and they 
have a linear spin wave dispersion relation with v s = \TplXs- The Euler-Lagrangian equations of motion gives the 
following dispersion relation for the n\ and n-i modes: 



,2 _ J 4 



\p 2 k A 



^fc 2 



2£ 

Ka 



for the AF state with (ni) = 0, (n 2 ) 7^ 0, 



for the SC state with (m) 7^ 0, (712) = 0, and 



±p<5fc 2 
^k 2 



UJ 2 = 



4(1 + sin 2 a/ cos 2 a)p 2 k 4 
2Sc °f a + (±5psm 2 a + -^)k 2 



(8.9) 



(8.10) 



*■ sin 2 a and (n|) 



for the mixed state with (n 2 ) — ,.„, ., 
lattice model results at the projected 50(5) symmetric point if we make the following identification 



s cos 2 a. These dispersion relations agree exactly with the 



2J 



Xs 



1 

2J 



2(4 J - A) 



i.ll) 



The effective Lagrangian can be easily used to discuss effects of 50(5) symmetry breaking. The simplest form of 
symmetry breaking is increasing the chemical potential beyond the critical value /i c , so that a pure SC state is 
realized. The chemical potential enters the effective Lagrangian through the gauge coupling in the time direction via 
the following substitution 



d t fii d t rii + 5pn 5 , d t n 5 d t n 5 — Spn 5 



1.12) 



where Sp = p — p c is the deviation of the chemical potential away from the critical value. In this case, the spin triplet 
excitations acquire a finite mass gap, with the following dispersion relation: 



■(k) = '— 

Xs 



pk 2 4(/i - p c ) 



Xs 



(8.13) 



and the mass gap increases with increasing doping in the SC state. 

We summarize the behavior of the collective modes obtained in the previous two sections in Fig. 



C0(k=0) 




C0~k 2 +A 2 



C0(k=0) 



00~k 



C0~k 



Pc P Pc P 

FIG. 5. Evolution of the collective mode spectra as a function of density in the projected SO(5) model. Fig. 5a shows the 
gap towards spin excitations. Charge excitations are gapless for the entire region of density, however, the dispersion relation 
changes from lo ~ k 2 to u> ~ k at p c , as indicated in Fig. 5b. 
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We see that while there are significant modifications of the collective mode spectra in the density regime < p < p c 
from the unprojected 5*0(5) symmetry, the spectra beyond p c is essentially identical to the behavior expected from 
the unprojected 5*0(5) symmetry. This should be expected from our general considerations about the Gutzwiller 
projection without much detailed calculations. We argued that the only effect of the Gutzwiller projection is to 
change the quantum commutation relation between the SC components of the superspin n\ and n$. However, for 
p > pc, the system is in a pure SC phase where these components acquire classical expectation values. In this case, 
the modification of the quantum commutation relation does not have any significant effect. This argument can also 
be illustrated by a simple picture of a chiral 50(5) sphere, as depicted in Fig. (^|). 



In this picture, the north and south poles represent the three AF directions, and the equatorial plane represent the 
SC directions. The sphere is perfectly 50(5) symmetric. However, the chemical potential along the pole direction 
acts like a fictitious magnetic field which restricts the sense of the rotation in the SC plane. Small oscillations of a 
vector pointing close to the north pole enclose the fictitious magnetic flux, and can only execute chiral rotations. This 
amounts to the projection of the particle-pair states at half-filling. On the other hand, small oscillations of a vector 
pointing anywhere along the equator does not enclose the fictitious magnetic flux, and their dynamics is therefore 
unaffected by the projection. Dynamics of a vector pointing anywhere between the north pole and the equator is also 
partially affected by the projection, but the symmetry of the static potential bears a unique signature. 



The main purpose of this paper is to introduce the concept of projected 5*0(5) models and discuss the properties 
of this model in connection with high T c superconductivity. The projected 50(5) model describes the low energy and 
long distance bosonic degrees of the freedom near the AF/SC transition. We showed that the Gutzwiller projection 
can be implemented analytically on every site in the 50(5) theory. In the presence of a infinite Mott-Hubbard gap, 
we show that static properties of the model can remain 50(5) symmetric, while the modification of the dynamics can 
be completely cast into a non-trivial commutation relation between the two SC components of the 50(5) superspin. 
Unlike the unprojected 50(5) models which can only have the full dynamic 50(5) symmetry at half-filling, the 
projected 50(5) model can have static 50(5) symmetry at a critical value of the chemical potential p, c and for a 
finite range of doping < p < p c . At p = 0, the system has a AF ground state and zero compressibility. In the 
intermediate regime < p < p c , the system has mixed AF and SC order and infinite compressibility. For p > p c , the 
system has a pure SC ground state, the SC order parameter rises to a maximal value before it decreases with doping. 
At the projected 50(5) symmetric point, we can understand precisely the evolution of the collective modes. On the 
AF side, we have two gapless spin wave modes and a gapless charge mode describing the gapless fluctuation from AF 




FIG. 6. Pictorial representation of a chiral 5*0(5) sphere. 



IX. CONCLUSION 
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to SC. In the intermediate density regime < p < p c , the physical properties of the spin waves remain unchanged, 
while the massive spin amplitude mode gradually decreases its energy and merges with the two spin wave modes at 
p = p c . The charge mode in the intermediate density regime is gapless, but has quadratic dispersion relation, which 
is a unique signature of the projected 50(5) symmetry. For p > p c , the charge mode is gapless with linear dispersion 
relation, and the tt triplet spin mode becomes massive, and gradually increases its energy with increasing doping. In 
this regime, the behavior of the collective modes are identical to the unprojected 50(5) model. 

This very simple model can form the basis to understand many novel and puzzling properties of the high T c 
superconductors in a unified framework. It points out a route from AF to SC through a gradual rotation of the 
superspin angle. At the projected 50(5) symmetry point the mean field energy is independent of the superspin 
angle, and therefore it offers an explanation of the absence of the chemical potential shift in the undcrdoped regime 
without global phase separation. It predicts a phase diagram which is qualitatively consistent with the observed phase 
diagram in the high T c materials. In the underdoped regime of the phase diagram, the system have large AF and SC 
fluctuations, and these fluctuations can be responsible for the pseudogap physics observed in these materials. 

There are many possible directions to carry out this line of research in the future. The most important issue is 
to understand the precise nature of the intermediate state in the regime < p < p c . Since the system has infinite 
compressibility in this regime, different small perturbation may select different ground states. Such perturbing effects 
might include quantum fluctuations and longer ranged interactions. In particular, we would like to investigate the 
possibility that these perturbations might lead to the formation of incommensurate order or stripes. 

In this work, we have discussed extensively the collective fluctuations in the long wave length limit. Due to the 
definitions of our effective lattice model, the k — ► limit corresponds to the k — > limit in the SC correlation functions 
and the k — > (tt, tt) limit of the AF spin correlation functions. Within the 5*0(5) theory, the tt resonance in the SC 
state is viewed as the 50(5) symmetry partner of the k — > Goldstone mode of the SC phase fluctuation. While the 
commensurate neutron resonance mode is observed in both YBCO and BISCO superconductors, all high T c systems 
also have incommensurate spin fluctuations. How can these features be explained within the current theoretical 
model? 

The fact that LSCO and YBCO have very different Fermi surface shapes and yet have similar incommensurate 
magnetic peaks strongly suggests that the incommensurate peaks are not sensitive to Fermi surface effects and should 
be explainable within an effective bosonic model. Let us recall that the collective mode of a superfluid boson system 
consists of linearly dispersing phonon branch and another roton branch with a minimum located at the inverse inter- 
particle spacing. So far, we have only studied the phonon branch of the charged bosons. By analogy, the roton branch 
should also exist, with a wave vector determined by the density of the charged bosons or doping. Within the 50(5) 
theory, while the commensurate neutron resonance can be viewed as the 50(5) partner of the SC phase mode, the 
incommensurate magnetic peaks can be viewed as the 50(5) partner of the roton minimum of the charged bosons. 
A detailed quantitative analysis of this picture will be carried out in the future. 

However, while the ground state in the doping range < p < p c may depend sensitively on small perturbation 
effects, at finite temperature, these perturbation effects should be small and the system should display more universal 
properties. We have shown that the projected 50(5) symmetry should be valid for the entire doping range < p < p c , 
and we shall quantitatively study the manifestation of this symmetry at finite temperature, and see if the projected 
50(5) symmetry can give a universal explanation of the pseudogap physics. 

Note Added: After completing this work, we received a very interesting paper by Coen van Duin in which he 
also observed the "remnant 50(5) behavior in the large U limit". 
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APPENDIX A: SLAVE BOSON RESULTS 



Alternatively, one can enforce the hard-core constraint (6.1) by introducing an additional "slave" boson for each 
lattice site. The presence of this boson (e(x)) indicates that the lattice site x is in the singlet state. The "less or 
equal" hard-core condition fl6.l|) is replaced with the equality constraint 



Q( x ) =^2 t l( x K(x) +t{(x)t h (x) + J(x)e(x) 



0, 



(Al) 



with q = I. Since in physical states one always has one and only one boson per lattice sites, destruction (creation) 
of a boson t a (i£) must always be accompanied by creation (destruction) of the empty boson (e). In this way, the 
physical operators for creating (destroying) a triplet (a = x,y,z) or a hole pair (a = h) acquire the form t\e (t a e'). 
The advantage of this method is that the constraint can be enforced exactly (at least in principle) by introducing 
an additional time-independent field A(x) at each lattice site, by adding to the Hamiltonian a term —\(x)Q(x) and 
by integrating over the X(x) on the imaginary axis. The whole Hamiltonian (2.11) thus takes the form (apart for a 
constant) 



H sb = - H X )Q( X ) + A sJ2 *l«(aO*a(aO + A c ^ t{(x)t h {x) 

x x,a x 

-J s /2 J2 {tUx)e(x) + h.c.)(tUx')e(x') + h.c.) 



(A2) 



<x,x'> ,a 



-J, 



./2 (t{(x)e(x)e f (x')t a (x')e(x') + h.c 



<x,x'> 



In practice, one starts with a mean-field approximation and expands the boson operators around their mean-field 



values as in Sec. VI. This expansion can be rigorously controlled by generalizing the constraint (Al) to large values 



of q, whereby one scales J c / S — > J c / s /s (cf. Refs. 1 33 34]). Physically, this corresponds to allow for a large number of 



bosons to be present at each site and thus to have a large value for the total spin, or more precisely for the SO (5) 
quantum number s, at each site. The mean-field result thus corresponds to the q — > oo limit, while the quadratic 
expansion corresponds to the first 1/q correction. 

At the mean-field level, the constraint is fulfilled exactly and indeed one obtains the same result and the same 
phase diagram as the variational ansatz (4.1) discussed in Sec. [v| By expanding the bosons quadratically around the 
mean-field one obtains the same modes as for the soft-core Hamiltonian with a similar dispersion. (Here, we restrict 
again to the SO(5)-symmetric case). Specifically, in the mixed phase one obtains two spin- wave modes with dispersion 



4 J 



k, 



one massive spin-amplitude mode 



and a quadratic mode 



u(k) 2 = (16J 2 - A 2 ) cos 2 a + 0(k 2 ) 



,{k) = 



4J + A fc4 



8 cos a 

In the pure superconducting phase, one has a 7r-triplet with dispersion 



and the SC Goldstone mode with dispersion 



w(k) 



, 16J 2 -A 2 j2 
- - k l 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 



These results coincide with the ones of the soft-constraint approximation Eqs. ( |6.6| ) to ( 6.1 C| ) in the limit of small 
4 J — A, i. e. , for low boson density. 
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